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Abstract. We search for nonlinear effects in nuclear giant resonances (GRs), in particular the isovector 
dipole and the isoscalar quadrupole modes. To that end, we employ a spectral analysis of time-dependent 
Hartree-Fock (TDHF) dynamics using Skyrme forces. Based on TDHF calculations over a wide range of 
excitation amplitudes, we explore the collectivity and degree of harmonic motion in these modes. Both GR 
modes turn out to be highly harmonic in heavy nuclei from A=100 on. There is no trace of a transition 
to irregular motion and multiple resonances are predicted. Slight anharmonicities are seen for light nuclei, 
particularly for 16 0. These are mainly caused by the spin-orbit splitting. 
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1 Introduction 

Nuclei display a variety of collective excitation modes, ro- 
tation, low-energy quadrupole or octupole vibrations, and 
giant resonances [TJ[2] • They share this feature with many 
other finite Fermion systems [3]. Nuclei provide a partic- 
ularly rich choice because they consist of two species of 
Fcrmions, protons as well as neutrons, and are fully self- 
interacting Fermionic systems without coupling to exter- 
nal fields (in contrast, e. g., to ions in metal clusters or the 
trap potential in quantum condensates). The low-energy 
modes in the MeV range are typically large amplitude 
modes with fission as an extreme example. Giant reso- 
nances (GRs), on the other hand, are considered to oscil- 
late with small amplitudes, although they are often very 
collective to the extent that they exhaust a huge fraction 
of the corresponding multipole sum rule. In that respect, 
they are the dominating modes in the nuclear excitation 
spectrum and have been intensively investigated over the 
past decades, for reviews see, e. g. [!J[5]. The present paper 
concentrates on GRs and tries to evaluate from a theoret- 
ical perspective the validity of the assumption of small 
amplitude motion for them. 

GRs are usually described by the well-known Random- 
Phase- Approximation (RPA) [TJ[2], which can be derived 
as the small-amplitude limit of time-dependent Hartree- 
Fock (TDHF) [2|S]. Most calculational schemes work out 
the linearization of TDHF explicitly, which leads to large 
reductions in computational expense for applications in 
restricted symmetries, see, e. g. 0[S]. For less symmet- 
ric systems, it is competitive and conceptually simpler to 
employ full TDHF and to run it practically in the linear 
regime by using only small excitations. That is the method 
of choice in molecules and clusters [9, 10] but has also been 



used successfully in deformed nuclei [TT]. In any case, the 
system is analyzed in terms of its spectrum of vibrational 
frequencies. That is, in fact, a classical concept. TDHF 
generally provides a classical description for the dynamics 
of a many-body system. Semi-classical re-quantization is 
needed to recover quantized excitation states from TDHF 
[r2 lll3lll4j . That is trivially and tacitly performed in RPA: 
small-amplitude motion is harmonic motion and the quan- 
tization of the harmonic oscillator identifies the vibra- 
tional frequency a; of a mode with its excitation energy as 
-E-exc = h^- The validity of that interpretation depends on 
the "harmonicity" of the mode. Truly harmonic motion is 
distinguished by the fact that the vibrational frequency is 
independent of the amplitude. It is the aim of this paper 
to check that assumption for nuclear GRs, in particular 
the isovector dipole GR and the isoscalar quadrupole GR. 

"Beyond RPA" effects have been studied before from 
different perspectives. Widely used is the hierarchical ex- 
pansion in n-particle-n-hole (nph) states. RPA emerges in 
lowest order at the level of lph (and reversely lhp) states 
2J. The next step is to include 2ph states, often called 
second RPA 15, 16, 17 . The coupling to the swarm of 2ph 
configurations has two effects: a shift of the resonance po- 
sition and a damping width, often called the collisional 
width [18]. The real shift of the peak and the width are 
typically of the same order of magnitude [19]. The mode 
remains nicely harmonic if shift and width are small, while 
larger values indicate nonlinear effects contributing to the 
observed peak energy. An advanced analysis of the ex- 
tremely dense 2ph spectra even allows looking for tran- 
sitions from regular to irregular motion in the nonlinear 
regime [20]. The extension to the 2ph space soon grows 
extremely expensive, however, and is thus often simplified 
by using rather small expansion spaces [161117] , schematic 
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interactions [21j , or pre-diagonalization by selecting a sub- 
set of coupled collective states [22ll23ll24] . 

An alternative access to anharmonicities is provided by 
TDHF. One merely has to watch the change of frequency 
of a mode with increasing amplitude [25j and in a later 
stage the onset of chaotic patterns [26] . Ideally, one should 
use periodic TDHF orbits to allow for a semi-classical 
quantization of TDHF states by computing the classical 
action along one closed orbit as a function of the amplitude 
[L2"irn?ll2"T] . A quantal eigenstate is encountered whenever 
the action crosses a value (n + \)h where n is some inte- 
ger. Strategies for finding branches of periodic trajectories 
have been successfully developed [28]. They are, however, 
very expensive to use for non-spherical modes and provide 
far too high a spectral resolution for our purposes. We 
are interested in a global test of harmonicity by watching 
the evolution of the average resonance peak with increas- 
ing amplitude. This is to some extent along the lines of 
[29] where a fluid-dynamical approach to TDHF was used 
for that purpose. Here, we are going to use fully fledged 
TDHF without any symmetry restrictions or constraints 
and invoke the systematic techniques of spectral analysis 
[9]fT0]fTT] carried forth into the nonlinear regime. This is 
an efficient method and allows to reveal the key features 
of a collective mode. That has been demonstrated, e.g., 
in a similar study of monopole resonances using the time- 
dependent relativistic mean-field model [30] , 

The results of that conceptually simple study will pro- 
vide insights in several respects. First, we obtain an esti- 
mate for the excitation energy of the quantized GR state. 
The difference between the excitation energy and the es- 
timate from the RPA frequency corresponds to the peak 
shift as obtained in the elaborate 2ph models. Second, we 
can preview the existence of double resonances for a given 
nucleus. Double resonances represent second excited states 
of harmonic oscillations and are conceivable for nicely col- 
lective and harmonic modes [j"7ll2"7ll3"T] . There is, in fact, 
experimental evidence for a double dipole resonance in 
heavy nuclei [32, 33J . Harmonicity when observed in the 
TDHF analysis at varying amplitude is a strong indicator 
for the existence of double, and higher, resonances. Third 
and finally, one may observe a transition to large dissipa- 
tion and chaotic motion when going more deeply into the 
nonlinear regime [26 , 30J . 



2 Formal framework 

We use the techniques of [II] . The TDHF equations are 
derived from the Skyrme energy functional including all 
terms, time-even as well as time-odd 34 . It has been 
worked out in recent studies of TDHF in the large am- 
plitude regime that a careful and consistent management 
of all terms, in particular including the time-odd contri- 
butions to the spin-orbit potential, is crucial to avoid arte- 
facts [35]. Wavefunctions and fields are represented on 
a three-dimensional Cartesian grid in coordinate space. 
The fast Fourier transformation is used to compute the 
action of the kinetic energy operator and to solve the 



Coulomb problem. The representation imposes no symme- 
try restriction whatsoever. This allows to explore the full 
phase space, particularly in a possibly nonlinear regime. 
A transition to chaotic motion would then be signified by 
dynamical symmetry breaking and widespread coupling 
to all conceivable modes of the system (see, e.g., [36]). 
The freedom is provided by our treatment but it has yet 
to be seen whether the system decides to go that chaotic 
path. The static mean-field equations are solved using an 
accelerated gradient iteration [37] . The ground state wave- 
function is excited by an instantaneous initial boost ap- 
plied in the same manner to each single-nucleon state, 

i.e. <pi gs ' ) — ► ip a = exp (ibD)ipa where D is the multi- 
pole operator associated with the mode of interest and b 
is the boost amplitude which will translate into the am- 
plitude of the mode with an excitation energy E* oc b 2 . 
The propagation of the TDHF equations is performed 
with exponential expansion of the time-evolution opera- 
tor exp(— ih m it) [38j . The expectation value of the multi- 
pole momentum D(t) — J2 a yPa(t)\D\'Pa('t)^ 1S recorded 

along the dynamical evolution. Finally, the signal D(t) is 
Fourier transformed into the frequency domain to obtain 
the spectral distribution of multipole strength. Some fil- 
tering is required to avoid artefacts in the spectra. We use 
windowing in the time domain with cos n (jrt/2Tj\ n ) and a 
large n to have a short effective time window with width 
of about 250 fm/c. This is chosen to smoothen spectral 
details and to concentrate the analysis on the average res- 
onance peak. Such a strong filtering allows also to use a 
conveniently sized numerical box [39j . We actually employ 
a 24x24x24 grid with spacing of 1 fm. 

There exists a great variety of Skyrme parametriza- 
tions which perform equally well concerning bulk prop- 
erties of the ground state but differ in details for more 
refined observables [40] . We find that all parametrizations 
which provide a reasonable description of GRs in heavy 
nuclei perform very similarly concerning (an-)harmonicity. 
Thus we confine the presentation here to the single parametriza- 
tion SLy6 from @T]. 

As test cases, we chose a selection of doubly magic nu- 
clei from 16 up to 208 Pb. Results for the isovector dipole 
GR are shown over the whole sample. The trends turn out 
to be very similar for the isoscalar quadrupole GR and the 
discussion is reduced to only two typical nuclei. 



3 Results and discussion 

Figure [T] shows the spectral distributions with varying ex- 
citation amplitudes for the test cases in our sample. Note 
the large width which was deliberately chosen to override 
the details of spectral fragmentation over the lph states in 
the vicinity of the GR (the analogue of Landau damping 
in bulk matter). That fragmentation width varies from 4 
MeV in 16 down to 2 MeV for 208 Pb. For spectral fil- 
tering, we use a fixed folding width of 1.6 MeV for all 
nuclei, which suffices to suppress the details and always 
deliver a unique and well-localized resonance. The plots 
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Fig. 1. Strength distributions for isovector dipole (left six panels) and isoscalar quadrupole (rightmost two panels) GR modes in 
six spherical nuclei. The modes were excited with an instantaneous boost at different boost energies as indicated. All calculations 
were done with the Skyrme parametrization SLy6. The threshold for neutron emission into the continuum, "n.cont." , is indicated 
in the plots. 



show a narrow frequency window about the resonance in 
order to better visualize the drift of the peaks with exci- 
tation amplitude. Other parts of the spectrum can safely 
be ignored, because there is not much strength outside 
the shown window. Within the plotting window, we see 
nicely concentrated resonance peaks in all cases. That fact 
alone already carries a message: the isovector dipole and 
isoscalar quadrupole resonances stay well inside the regime 
of regular motion for all amplitudes considered here, and 
probably far beyond that point. That agrees with the find- 
ings of [29] • Note that the isovector monopole mode was 
found to be closer to the chaotic regime [3U] and that the 
highly nonlinear low-energy modes also very quickly turn 
chaotic [26] , 

The feature of interest for our studies are the trends 
with increasing amplitude. Reading the left three panels 
in figure [1] from bottom up, it becomes obvious that small 
nuclei show some changes while the spectra of heavy nuclei 
stay extremely robust. The middle panel checks the effect 
along an isotopic chain, here the Sn isotopes. Note that 



this represents not only isotopic trends but also dramatic 
differences in overall binding with 100 Sn being at the edge 
of stability, 132 Sn on the neutron rich side not much better 
bound, and 120 Sn well inside the valley of stability. There 
is very little difference in spectral pattern in view of these 
dramatic variations in isospin and stability. System size 
seems to be the key feature which determines the robust- 
ness of the resonance. At second glance, one finds that the 
two unstable nuclei (upper middle and lower middle pan- 
els) indicate indeed slightly larger drift of the resonance 
peaks than the stable 120 Sn. But this remains a faint ef- 
fect. As a rule of thumb, we find that the dipole resonance 
is a harmonic mode in nuclei with A > 100. 

The right panels in figure [1] show two samples for the 
isoscalar quadrupole mode. The changes are still small but 
generally stronger than for the dipole modes in the same 
nuclei. Note that the stronger shifts develop towards the 
lower end of the spectrum. This is caused by the much 
more nonlinear low-lying quadrupole mode [26] which, al- 
though having comparatively small strength, affects the 



4 



P.-G. Reinhard et al.: Nuclear Giant Resonances and Linear Response 





13.2 




14.8 


> 

0) 

S 


13 




14.6 


>. 
o 


12.8 




14.4 


requei 


12.6 




14.2 


eak f 


12.4 




14 


Q. 


12.2 


208 Pb dipole 


13.8 



10 



20 



25 





18.4 


> 


18.2 


m 




S 


18 


C 


O 
3 




cr 


17.8 










o 


17.6 


Q- 






17.4 




20 


> 


19.8 


0) 






19.6 


>. 
o 


19.4 


c 




<u 

3 


19.2 


cr 




0) 


19 


TO 


18.8 


O 




Q. 


18.6 




18.4 







40 Ca dipole I 





5 10 15 20 

excitation energy E [MeV] 



15 
14.8 
14.6 
14.4 
14.2 

14 







132 Sn dipole 




5 10 15 20 



120 



Sn dipole 



14.2 
14 
13.8 
13.6 
13.4 



8 10 12 14 16 




1 20 

Sn quadrut 


)ole 






5 10 15 20 25 







21.4 










21.2 










21 










20.8 






100 Sn dipole 




20.6 


16 quadrupole j 




5 10 15 


20 




D 5 10 15 20 


25 



excitation energy E [MeV] 



excitation energy E [MeV] 



Fig. 2. Trend of resonance peak frequency with excitation energy for isovector dipole and isoscalar quadrupole resonances 
in various spherical nuclei as indicated. All calculations were done with the Skyrme parametrization SLy6. The diagonal of 
"peak-frequency = excitation-energy", i.e. i? GXC = ftui, is indicated by a faint straight line. In case of 208 Pb, we also show a 
second steep line which stands for the the second resonance mode with E exc = 2hw. In case of ie O, results from a calculation 
without spin-orbit force are indicated by open circles. 



lower side of the resonance peak. Note also that the reso- 
nance as such is much more robust and the drift of peak 
position is not much larger than in the dipole case. Thus 
we can generalize the conclusion from above and call also 
the isoscalar quadrupole resonance a harmonic mode for 
medium-heavy and heavy nuclei. 

Figure summarizes the trends of peak positions with 
excitation energy. It confirms what we could read off (less 
quantitatively) from the previous figure. The dipole reso- 
nance is harmonic for A > 100. The quadrupole still shows 
slightly more drift in 120 Sn, but also that can be consid- 
ered as small and is further reduced for larger nuclei. 

The figure also shows faint steep lines which indicate 
the simple-minded oscillator quantization -Eoxc = The 
point where both lines cross each other indicates the am- 
plitude of the oscillation which contains precisely one os- 
cillator quantum. It is practically the oscillation frequency 
in the linear regime (extremely small amplitudes) for the 
harmonic modes in A > 100. There are sizeable correc- 
tions for the lighter nuclei. The actual resonance peak 
is 0.6 MeV below the RPA frequency for 16 and 0.25 
MeV below for 40 Ca. For the particularly harmonic case 
of 208 Pb, we also analyze the occurrence of the double res- 
onance excitation, see the second steep line in the upper 
left panel of figure O The mode remains safely harmonic 
up to the double resonance regime which confirms once 



more the existence of a double resonance without any fre- 
quency shifts |32j . 

It is known that the spin-orbit force has a tendency 
to degrade collectivity and enhance dissipation [32] • We 
have tested that for our worst case of 16 O. The lower left 
panel of figure [2] shows also results from TDHF calcula- 
tions without spin-orbit force. We see indeed more har- 
monicity than for the full force. This confirms again that 
the spin-orbit force to some extent counteracts collectivity 
in GRs. The mechanism is that the spin-orbit splitting en- 
hances spectral fragmentation which, in turn, reduces the 
coherence in the superposition of the many lph modes 
to one dominant collective mode [TJ[2] . Without spin-orbit 
force, GR would be highly harmonic modes all over the 
periodic table including light nuclei. 



4 Conclusion 

We have analyzed the isovector dipole and isoscalar quadrupole 
giant resonance (GR) in nuclei of different size using fully 
fledged time-dependent Hartree-Fock (TDHF) dynamics 
in connection with Skyrme forces. The aim was to check 
the commonly accepted view that GRs are collective modes 
corresponding to rather harmonic oscillations of the nu- 
clear many-body system. To that end, we analyzed the 
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spectral distributions of the corresponding multipole mo- 
ments for widely varying amplitudes of the oscillations 
and for a broad range of nuclei between the rather small 
16 and the heavy 208 Pb. We find in all cases that the 
GRs remain in the regime of regular motion confirming 
from that aspect the robust collectivity of these modes. 
"Harmonicity" is more critical. We see it extremely well 
fulfilled for heavy nuclei, at least from A = 100 on. This, 
in turn, suggests the existence and robustness of highly 
collective multiple-resonance excitations which have been 
experimentally assessed already for the double resonance 
in heavy nuclei. Light nuclei show first glimpses of non- 
linear effects to the extent that the resonance-peak po- 
sition shows visible drift with increasing amplitude. For 
16 0, e.g., we observe a drift of 0.6 MeV from the strictly 
linear regime (small amplitudes) to the realistic amplitude 
where the oscillator quantization recipe -E OS c = fru is sat- 
isfied, ft is still a small, but non-negligible, effect. The 
isoscalar quadrupole mode is similarly robust concerning 
the peak drift. The patterns of the spectral distribution, 
however, are a bit more volatile at the lower tail, which 
is due to the perturbing influence of the highly nonlinear 
low-energy collective surface quadrupole mode. 
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